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We obtain the equations of motion for cosmic strings in extensions of the 3+1 FRW model with 
extra dimensions. From these we derive a generalisation of the Velocity-dependent One-Scale (VOS) 
model for cosmic string network evolution which we apply, first, to a higher-dimensional isotropic 
D + 1 FRW model and, second, to a 3-1-1 FRW model with static flat extra dimensions. In the former 
case the string network does not achieve a scaling regime because of the diminishing rate of string 
intersections {D > 3), but this can be avoided in the latter case by considering compact, small extra 
dimensions, for which there is a reduced but still appreciable string intercommuting probability. 
We note that the velocity components lying in the three expanding dimensions are Hubble-damped, 
whereas those in the static extra dimensions are only very weakly damped. This leads to the 
pathological possibility, in principle, that string motion in the three infinite dimensions can come 
to a halt preventing the strings from intersecting, with the result that scaling is not achieved and 
the strings irreversibly dominate the early universe. We note criteria by which this can be avoided, 
notably if the spatial structure of the network becomes essentially three-dimensional, as is expected 
for string networks produced in brane inflation. Applying our model to a brane inflation setting, 
we find scaling solutions in which the effective 3D string motion does not necessarily stop, but it 
is slowed down because of the excitations trapped in the extra dimensions. These effects are likely 
to influence cosmic string network evolution for a long period after formation and we discuss their 
more general implications. 
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I. INTRODUCTION 

There has been a recent resurgence of interest in cosmic strings both for theoretical and observational reasons. Of 
particular interest is the genericpossibility that cosmic strings can be produced at the end of an inflationary phase 
in models of brane inflation 0, Q| . The evolution of the network of the strings created in these models can be very 
different from the standard field theory case, thus providing a potential observational window on superstring physics 
[^^1^. On the observational side, there are perennial cases of astrophysical phenomena for which cosmic strings have 
been invoked as an explanation in the absence of some more orthodox mechanism; recent examples are two peculiar 
gravitational lensing events 0,|3iSI3 but which need further independent follow-up. However, expected improvements 
in observational data, particularly from high resolution CMB experiments and gravitational wave detectors, present 
us with the very real prospect of detecting or constraining cosmic strings over a wide range of predicted energy 
scales (see, for exam ple , 0)^3)- Further theoretical motivation has come from a recent phenomenological study of 
SUSY GUT models [Ul, which again found generic cosmic string production (for all cases which solve the monopole 
problem). 

It is the recent work on cosmic strings appearing in higher dimensional theories, such as brane inflation, which 
primarily motivates our present study. For spacetime dimension greater than four, strings no longer generically 
collide, so that loop production will be highly suppressed. Loops radiate away energy from the long string network, so 
this suppression will result in a much higher density of cosmic strings than in the usual 3 -f 1 dimensional case. Jones, 
Stoica and Tye p| have estimated this enhancement by using a three-dimensional one-scale model and introducing an 
intercommuting probability P < 1 to account for the fact that strings generically miss each other due to the presence 
of extra dimensions. They suggest that the enhancement on the energy density of the string network is of order P~^, 
which can be orders of magnitude different than the usual case. 

This approach however does not take into account string velocities in the extra dimensions. In general, cosmic 
strings are subject to the constraint that the average velocity squared of string segments must be less than 1/2. Thus 
the fact that strings are moving in the extra dimensions will slow down their apparent three-dimensional motion. One 
might naively expect that velocities in the infinite dimensions will be redshifted by the expansion, while velocities 
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in the compact dimensions will not if these dimensions are static. Hence there is the cosmologically dangerous 
possibility that velocities in the extra dimensions will accumulate and dominate while string motion in the three 
infinite dimensions will come to a halt. This can occur before SUSY breaking, when the fields corresponding to the 
string position in the extra dimensions are expected to become massive. However, if 3D string motion stops for long 
enough in the early universe, then the strings would irreversibly dominate the energy density of the universe making 
a subsequent hot big bang model impossible. Thus we need to go beyond the simple analysis presented to date in 
order to gain a better quantitative understanding of cosmic string evolution in higher dimensional spacetimes, notably 
taking into account the important role of velocities in the extra dimensions. 

String evolution in three spatial dimensions has been studied by various authors. Kibble ^3] described string 
networks by a single lengthscale, the 'correlation length', and showed that it evolves towards a scaling solution in 
which it stays constant with respect to the horizon size. Bennett [T^ IT3 | later modified this 'one-scale' model with 
similar conclusions, subject to a condition on the efficiencyof small loop production. The existence and stability of 
this scaling solution was verified by numerical simulations [l^ [l^ • These studies also revealed new physics at smaller 
scales, in particular the accumulation of significant small-scale structure on strings, which results in loop production 
at much smaller scales than initially thought. 

To try to incorporate small-scale structure in analytic models, a number of different approaches have been attempted. 
These include a 'kink-counting' model [TtI IT^. a functional approach 19j, a 'three-scale' model 20] and a 'wiggly' 
model |2lj |. Including small-scale structure in analytical models comes with the cost of introducing several extra 
parameters, which need to be fixed by simulations. 

However, the large-scale properties of string networks can be quantitatively described by the Velocity-dependent 
One-Scale (VOS) model 22, 23,'23|, which does not suffer from this problem. By introducing a variable rms string 
velocity, the VOS model extends its validity from the friction dominated regime at early times, through the matter- 
radiation transition to A-domination at late times , thus describing the complete cosmological history of string net- 
works. Though it does not directly model small-scale structure, it provides a 'thermodynamic' large-scale description 
of cosmic string evolution, which agrees remarkably well with high resolution numerical simulations. 

The purpose of this paper is to extend the VOS model to spacetimes of higher dimension. Although strong 
motivation is provided by brane inflation, where the extra dimensions are small and stabilised, we intend to keep the 
discussion as general as possible to include time-varying extra dimensions, as for example in a (Z) + l)-dimcnsional 
FRW universe. 

The paper is organised as follows. In section^lwe discuss cosmic string dynamics in a FRW spacetime with isotropic 
and flat (but possibly expanding) extra dimensions. Starting from the Nambu-Goto action in this spacetime we derive 
the equations of motion as well as an expression for the energy of a cosmic string network. In section IhH we derive 
the averaged equations describing the evolution of cosmic strings in that {D + l)-dimensional spacetime. After we 
briefly review the 3-1-1 VOS model, we extend it to Z? -I- 1 dimensions and comment on some qualitative features of 
solutions. In section Hvl we consider possible application of this extra-dimensional VOS (EDVOS) model to the case 
of brane inflation and discuss the dependence of the results on various parameters as for example the intercommuting 
probability of strings. We conclude in section^ Finally, there is an Appendix where we derive approximate formulae 
for the momentum parameters of the (D + l)-dimensional VOS model. 



II. COSMIC STRING DYNAMICS IN HIGHER DIMENSIONS 



We consider a cosmic string propagating in a I? + 1 dimensional spacetime with metric g^i, (/i, ^' = 0, 1, 2, . . . , D). In 
the limit that its thickness is much smaller than its radius of curvature, the string can be regarded as a one-dimensional 
object with a world history described by a two-dimensional spacetime surface, the string worldsheet 

xf=x''{C), a = 0,1. (1) 
The dynamics is given by the Nambu-Goto action 




where fj, is the string tension and 7 is the determinant of 7^^ — g^ydax'^df^x'^ , the puUback metric on the worldsheet. 
The equations of motion for the fields x^ obtained from this action are given by 

+ T''^^-i°'Pda,x''dpx^ = (3) 

where F^^^^^ is the {D -\- l)-dimensional Christoffel symbol 

r^A = ^5'^" (SAffK. + a,g,A - (4) 



and V^x^ the covariant Laplacian of the worldsheet fields x^^ given by 



(5) 



By varying the action with respect to the background metric we obtain a spacetime energy-momentum tensor 



(6) 



We wish to study the evolution of cosmic strings in a FRW universe with D — 3 extra dimensions. For simplicity 
we choose the following metric, allowing toroidal compactification of the extra dimensions 



ds^ = N{tfdt^ - a{tfdx^ - b{tfd\^ 



(7) 



where t = x^, x = a;' with i = 1,2,3 and \ = x^ with t = 4,5,..., D. The lapse function N{t) allows us to switch from 
cosmic {N{t) = 1) to conformal time by simply setting N{t) =a{t). For stabilised compact extra dimensions, b{t) is 
set to a constant and the coordinates 1 are periodically identified. Alternatively, the scalefactor b{t) of the '€-space' 
can be set to a{t) for a generalised {D + l)-dimensional FRW universe or even more generally be allowed to evolve 
independently of a{t). 

The action |(5Jl is invariant under worldsheet reparametrisations, which we can use to choose a gauge. In flat space 
either the 'light-cone' = t, = z — t) or the 'conformal' gauge (i ■ x' — 0, x'^+x' =0) is usually chosen but for 
cosmological backgrounds it is convenient to work in the gauge 



C° = i, x-x' ^0 



(8) 



where dots and primes denote derivatives with respect to the timelike and spacelike worldsheet coordinates and 
respectively. By choosing this gauge we identify spacetime and worldsheet times, while we impose that the vector x 
is perpendicular to the string tangent, thus representing the physically observable velocity. 
In this gauge the equations of motion ^ in the spacetime ^ are 
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where e is a scalar, the energy per unit coordinate length, defined by 
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The energy-momentum tensor (jSJ for the metric Q in the same gauge becomes 

= Nu^Id-s i^i'i" - e-'x'^^n <5(^) ( x - x(C, t),l- 1(C, t) ) 

where we have defined ^ = C^. 

The energy of the cosmic string can be defined from T'^" as 



E= f Vhni,n^T^"'d^iid°-H 

J i^const 
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where = {N, 0) is the normal covector to the spacelike D-dimensional surface t = const and h is the determinant 
of the metric on that surface related to the volume element by 
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I'^dx''^ A dx>'^ A... dx^" = Vhd^xd^-^l . 



(15) 



With equation H13|l the energy becomes 



E{t) = N{t)nJ ed(. (16) 

Note that when 1 = and N{t) = a{t), equations (|9I11|I and H16|l reduce to the usual equations of motion and energy 
definition for a string in an FRW universe written in conformal time ,25j. These have been used to study cosmic 
string evolution in an expanding (3+l)-dimensional universe. In the next section we discuss cosmic string evolution 
in higher dimensions, based on equations H9I11|) and Hlt)|) . 



III. COSMIC STRING EVOLUTION 



In this section we discuss the cosmological evolution of strings. After reviewing some basic methods and results in 
three dimensions, we derive the equations describing the evolution of cosmic strings in a higher dimensional spacetime 
with metric 0, giving particular attention to the case of an isotropic, (Z) + l)-dimensional FRW universe {b{t)—a{t)), 
and that of stabilised extra dimensions b{t) — const. Starting from equations (|9I1H) and (|16|) we write down a 
higher dimensional extension of the Velocity-dependent One-Scale (VOS) model j2^, which has been successfully 
used to provide an analytic 'thermodynamical' description of the basic properties of evolving cosmic string networks. 
Application to cosmologically interesting cases is also discussed. 



A. Basics: Evolution in 3+1 Dimensions 



Monte-Carlo simulations of cosmic string formation after symmetry breaking phase transitions suggest that to a 
good approximation the strings have the shapes of random walks at the time of formation. Such 'Brownian' strings 
can be described by a characteristic length L, which determines both the typical radius of curvature of strings and 
the typical distance between nearby string segments in the network. On average there is a string segment of length 
L in each volume and thus the density of the cosmic string network at formation is 

P-JJ-Jl- (17) 

A heuristic picture of the evolution of the density of the string network can be obtained as follows. Assuming that 
the strings are simply stretched by the cosmological expansion we have p ~ a(t)~^. This decays slower than both the 
matter and radiation energy densities and so such non-interacting strings will soon dominate the universe. 

This picture changes when the effect of interactions is taken into account. As the network evolves, the strings 
collide or curl back on themselves creating small loops, which oscillate and radiatively decay. Via these interactions 
energy is lost from the network so that string domination can be avoided. Each string segment travels on average a 
distance L before encountering another nearby segment in a volume L^. Assuming relativistic motion and that the 
produced loops have an average size L, the corresponding energy loss is given by pioops ~ L^'^jiL. The energy loss 
rate equation becomes 

p«-2^p-f (18) 

Cosmic string networks are known to evolve towards a 'scaling' regime in which the characteristic length L stays 
constant relative to the the horizon du ~ t To see this we set L = ^{t)t and substitute l(T7j> into l(TH|> to obtain 

^ - ^ ^2{P-l) + -\ . (19) 



7 2t V ' '7 

The parameter (i is related to the scalefactor a[t) by a[t) oc and is equal to 1/2 and 2/3 in the radiation and matter 
eras respectively. Equation (|19|l has a scaling solution 

7=[2(l-/3)]-i (20) 

(which depends on cosmology through the expansion exponent (3) demonstrating that the characteristic length scales 
at a value L ^ t. If we start with a high density of strings, intercommuting will produce loops reducing the energy of 
the network, whereas if the initial density is low then there will not be enough intercommuting and 7 will decrease. 



Given enough time, the two competing effects of stretching and fragmentation wiU always reach a steady-state and 
the scahng regime will be approached. 

Equation (jl9|l was derived on physical grounds and it only captures the basic processes involved in string evolution, 
namely the stretching and intercommuting of strings. It does not take into account other effects like the redshifting 
of string velocities due to Hubble expansion. However we can derive a more accurate evolution equation for the string 
energy density based on the Nambu-Goto action. By differentiating the energy ()16|1 with respect to time for the case 
of a 3+1 FRW universe and setting E oc pa"^ we find 



2- + 2N-^aa{±^)] p (21) 



where we have defined 

We also introduce a phenomenological term to account for energy loss through loop production 

Ploops — ^ 

where c is the loop production parameter, related to the integral of an appropriate loop production function over all 
relevant loop sizes (see equation (|37|) and the discussion above it), and v the average velocity of intercommuting string 
segments. 

Using (|17(l and setting L = j{t) t as before, we obtain the following equation, expressed in physical time t 

7"^$-:if2/3(l + .^)-2+^V (24) 



dt 2t \ ^ ' 7 
Here we defined w, the average velocity of string segments, by 



c?T / /et^C 

with T the conformal time {N(t) — a{T)). 

Equation (|24() is of the same form as (|19() but has an extra correction term (3v'^ accounting for redshifting of velocities 
due to cosmological expansion. It also includes the parameter c, the value of which can be extracted from numerical 
simulations and it is of order unity [23, • 

To solve (I24|l we also need an evolution equation for v. This can be obtained by differentiating (|25|l with respect 
to T and using the three-dimensional version of the equation of motion for x HlOf) . The result (expressed in physical 
time t) is ^ 

= (1 - .^) f 4 - 2^^^') (26) 



dt ' ' \R 

where H is the Hubble parameter a/ a = (3t~\ R the average radius of curvature of strings in the network and k the 
momentum parameter (first introduced in |22|) defined by 

kv{\ — v"^) , . ^, , . 

^ ^ = (i . u(l - x)) (27) 
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u being the curvature vector defined by 



d^x 1 . 



with ds — Vx^ dC,. For a Brownian network (and within the VOS assumptions) the average radius of curvature R is 
equal to the correlation length L = jt. 

Equation H2t)|) has a clear physical meaning: velocities of string segments are produced by string curvature and 
damped by the cosmological expansion. Together with H24() it constitutes the Velocity-dependent One-Scale (VOS) 



model, which has been demonstrated to be in very good agreement with numerical simulations 26] . It has the scaling 
solution 
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(29) 



in terms of the expansion exponent /3, the loop production parameter c and the momentum parameter k. 

The momentum parameter is a measure of the angle between the curvature vector and the velocity of string segments 
and thus is related to the smoothness of the strings. Slowly moving strings are smooth so the velocity is more or less 
parallel to the curvature vector, resulting in a value of k of order unity. As v increases towards relativistic values the 
accumulation of small-scale structure renders the strings wiggly. Velocities become uncorrelated to curvature and k 
decreases. In particular it can be shown analytically that for flat space, where = 1/2, the momentum parameter 
vanishes for a wide range of known solutions [22, 123 ■ 

An accurate ansatz for the momentum parameter k has been proposed in [23 | 



k — k{v) 



2V2 1 - 8v^ 

TT 1 + 8t)6 



(30) 



satisfying k{l/V2) = 0. 

Note that the fact that v — l/\/2 in flat spacetime, can be shown analytically for closed loops only, but for long 
strings it is observed in numerical simulations ^25]. For expanding or contracting spacetimes, v is less or greater than 
l/\/2 respectively. Hence for an expanding universe, string velocities are subject to the constraint 

v'<l. (31) 
This fact will be important for our discussion of extra dimensions. 

B. Evolution in Higher Dimensions 

We now proceed to derive macroscopic evolution equations for string networks in higher dimensions, based on 
equations (|9I11|I and H16|l. derived in sectional The result is the extra-dimensional velocity-dependent one-scale 
(EDVOS) model. 



1. Evolution Equation for 7 



We start by differentiating the energy equation with respect to time and using the equation of motion for 
e. This gives 
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where we have defined the average of a function / along the string as 
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We then use E cx pa^b^ and definition H12|l for e to obtain an evolution equation for the energy density of the string 
network 
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We would normally proceed by assuming that the string network at formation is Brownian with a characteristic length 
L and writing 



PlL jjL 



(35) 



Unfortunately, since the metric we are considering is not in general isotropic, we do not expect this Brownian structure 
to be preserved by the evolution. The amount by which string segments are stretched will depend on their orientation 
and as time passes, the random walk shape of strings will be distorted by the anisotropic expansion. However, there 
are two interesting special cases in which the strings can remain Brownian. First, the case of isotropic expansion 
{a{t) =b{t)), corresponding to a generalised {D + l)-dimensional FRW universe, and second, a situation where the 
formation of the string network is localised on an isotropic slice. The latter is the case, for example, in brane inflation 
p7t l28t l29l Isd l3ll | where string formation takes place essentially on a brane and there is an effective three-dimensional 
description of the evolution. We will consider both cases below. 

(i) Isotropic case: Starting with the isotropic case, we set b{t) = a(t) in equation (|34|l to obtain an evolution equation 
for the energy density of a non-interacting string network in a (D + l)-dimensional FRW universe 



(D- 1) -f 27V"^a^(x^) -f 27V"^a^(r) p. (36) 



As in the usual VOS model we can treat string interactions by introducing a phenomenological loop production term. 
In the three-dimensional case a string segment of size L, the correlation length, is expected to travel a distance L 
before encountering another segment and interacting with it in a volume L^^. If the probability of such an interaction 
producing a loop of length between £ and -t- c?^ is given by the scale invariant function f{£/L) then the energy loss 
due to loop production is j2j| 



HV de vp f°" dl dvp 

Pioops^-;^/ -^mL) = -j; -mL)=-—. (37) 



When the number of spatial dimensions D > 3, the string segments will generically not interact after moving a 
distance L. In particular, if the strings have a thickness or capture radius 6 (quantifying how close they need to be 
in order to interact), then the probability of interacting after moving a distance L (in time 6t) in a volume is 



vSt 1 /S^ 



(38) 



Note that this suppression is essentially the intercommuting probability P which we shall further discuss later, though 
here it is time dependent. The loop production term is then 

/ c \ -D — 3 /.oo 7/7 ~ / r \ 

pv I \ / j./„ /r^ I " \ ~ 



= -I5(lJ i. T^Wi) = -T^Uj ^-f- 

where we have defined 



CZ5 = 5 j . (40) 

Including this loop production term in equation II36I) and setting as in the three-dimensional case L = ^{t)t, where t 
is the physical time {N{t) — 1) we obtain 



-1^7 1 
7 -rr^ 



dt {D-l)t 



P{{D-l) + 2v^) + ^-^-{D-l) 



(41) 



Here = v^"^ -t- where Vx and vi, the rms peculiar velocities of string segments in the x and the £ directions 
respectively, are defined by 



since b(t)=a{t) in the isotropic case. 

Note that since cd is proportional to L^*^^^^^ — {■-ft)^'^^^^\ the loop production term vco/^ is explicitly time 
dependent and no scaling solution 7 = const exists. In particular, as time increases the loop production term becomes 
smaller and smaller, reflecting the fact that strings cannot find each other and interact in more than three spatial 
dimensions. String evolution in this isotropic D + 1-dimensional case will be discussed further in section UlIB 41 

(ii) Anisotropic case: We now consider the case were b{t)^a(t) in lO but the string network is produced on a FRW 
slice (the x-space) of the spacetime. If the extra dimensions are larger than the correlation length and the strings 



are allowed to move in the bulk after formation, then the anisotropic expansion will distort their Brownian structure 
and they will be more difficult to characterise analytically. We leave this case for subsequent discussion (see the final 
section). The situation is more straightforward if the extra dimensions are compactified at a size smaller than the 
correlation length which we shall assume from now on. We can then have an effective three-dimensional description 
in which the strings will remain Brownian due to isotropy in the three-dimensional FRW space. The presence of the 
extra dimensions will reduce the intcrcommuting probability in the effective three-dimensional description, but we 
shall assume it is non-zero. In principle, we can also have time- varying extra dimensions, as long as their size remains 
much smaller than the correlation length at all times. However, for simplicity we only consider static extra dimensions 
and set b(t) = l, assuming stabilisation by some unspecified mechanism. 

For such cases we can substitute (jl7|> in equation H34|l to obtain an effective three-dimensional evolution equation 
for the energy density, which however takes into account velocities in all D spatial dimensions. As before we also 
introduce a loop production term — and set L = 7(t) t where t is the physical time {N{t) = 1). We thus obtain 
an effective evolution equation for j{t) 

7"' ^ - ^ |/3 [(2 + zi;|) + (2 - wj) + (l - w',) v,'] - 2 + ^| . (43) 



In deriving this expression, we have assumed 
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and we have defined an orientation parameter 



which quantifies the degree to which the strings lie in the extra £ dimensions. 

Note the differences between (I43II and the analogous equation for a 3+1 FRW universe H24|l . Comparing 1431) to H24|) 
we see that the coefficients of dilution and a;-velocity redshift are modified by w'j , while there is an additional term 
accounting for the redshifting of velocities in the static ^-space. We interpret this term as arising from the conservation 
of momentum: an expanding string segment which is transverse to the i'-space and moves in the i directions is slowed 
down because its effective mass increases as the string is stretched in the x-direction. 

Another important difference with the three-dimensional case is that the loop production parameter c is now 
suppressed by a factor of order {5 / Ri)''^~^\ where Ri is the compactification scale. If Ri was greater than the 
correlation length L, then the suppression factor would instead be {5 / LY^~^\ which is time-dependent. Thus there 
would be no scaling solution initially, as in the (D-|- l)-dimensional FRW case we discussed above (see later discussion 
for the asymptotic evolution in this case). Here however, we assumed the compactification scale is much less than 
the correlation length and so the suppression factor can be treated as a constant intercommuting probability. For a 
string-theoretic calculation of this probability see Ref. |^. 



2. Evolution Equations for Vx and ve 



As in the three-dimensional case we also need to know how the average velocities of string segments evolve in time. 
We begin by considering the case of static extra dimensions b(t) = 1. Working in conformal time r {N{t) = a(T)), 

we start by differentiating Vx"^ = with respect to r and using the equation of motion H1Q() for the fields x. We 

obtain 
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where u is the curvature vector defined by 
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(47) 



with ds ~ y (x'^ + l'^/a^)(iC £ind R the radius of curvature, equal to the correlation length L = for a VOS model 
Brownian network. 

The first three terms in equation H46(l (first line) can be neglected because their only effect is to slightly modify the 
coefficients of all the remaining terms in the expression (but for the last). Numerical confirmation of the smallness of 
such terms (in particular the first term in the case of conventional FRW) was presented in [2^ . The remaining terms 
(again except the last one) are of the same form as in the (3+l)-dimensional FRW case but now they also include 

terms dependent on the velocity in the ^-space and an extra term which gives corrections of order ' 



The last term arises from the artificial splitting of the Z?-dimensional velocity into an x and an H part. It is 
absent in the usual (3+l)-dimensional case because of the gauge condition x • x' = 0. Here the gauge condition is 
X • x' — 1 • I'/a^ = and thus x • x' does not identically vanish. However, we expect its value along the string to change 
sign randomly with no large-scale correlations. Hence, averaged over the whole string network, x • x' is expected to 
vanish so that we can neglect the last term of H46|l . This was tested numerically for a three-dimensional FRW model, 
splitting V into a xy- and a z-part and neglecting the corresponding terms in the xy- and z-velocity equations. The 
evolution of the system using these equations was found numerically to be very close to the corresponding evolution 
using the full three-dimensional velocity equation, with the scaling values of v and 7 agreeing to three significant 
figures. 

Again using the approximation H44|l and switching to physical time i, we can write (|46|) in a much more elegant 
and useful form, 

^^dv^^ _ ^2) _ ^2 „ ^^^2(1 _ ^2) _ ^^^2^^2 (48) 

where wi is given in (|45|l and kx is defined by 
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Similarly the evolution equation for V( reads 

~dt ~ ^ 
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with ki defined in analogy to k^ 
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We see that, as may have been anticipated from (|34|l . the Vx damping term is very different from that of vg. To 
demonstrate this we neglect w'j corrections and note that Hvx^ in H48(l comes with a factor of (2 — 2vx^ — vg') but 
Hvi' in (|5()|l has a factor of only (1 — — vi') (which can cancel almost to zero). This result will be important for 
the discussion in the next section. 

We can also write down an evolution equation for the total velocity u, by differentiating = ((x^, 1^/a^)). The 
result is simply the sum of (|^ and (|^ with 

kxVx + kivi = kv (52) 

where k has a similar definition to l|49|l and (|51|l involving the dot product of (x, 1/a^) with the curvature vector u. 
Equation H52|) follows from the linearity of the dot product and the integral. 

The momentum parameter k, as in the (3-|-l)-dimensional case measures the angle between the curvature vector 
and the velocity of string segments, thus providing a measure of the smoothness of the strings. For smooth strings 
(the non-relativistic limit) the velocity of string segments is expected to be more or less parallel to the curvature 
vector u corresponding to k of order unity. For relativistic velocities however, small-scale structure accumulates and 
the strings become wiggly. The parameter k is then expected to approach zero. 

Similarly, kx and ki provide a measure of the angle between the curvature vector and the x- velocity or the velocity 
respectively. Therefore they encode two effects: the wiggliness of the strings and the extent to which the curvature 
vector u lies in the x or ^ subspace. 

Following [231 we split v = (x, i/a) into a 'curvature' component Vc produced during the last correlation time and a 
'left-over' component Vp, coming from previous accelerations. We can then derive the following approximate formulae 



for k,kx and kg, valid in the relativistic regime (see Appendix 
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(55) 



Note that for fiat spacetime, the condition v'^ — \/2 still holds and so fc(f/\/2) =0. For an expanding universe the 
velocities of string segments are subject to the constraint 

v^=vl+vj<]^. (56) 

Finally, we note that the velocity evolution equation in the case of an isotropic [D + l)-dimensional FRW universe 
{h{t) = a{t)) can be similarly found 

| = (l-„=,(i-2H..). ,57) 

Note that this equation does not depend on wi. However, in the case h{t) = 1 (eqns (|48|l . (|50|l 'l there is an explicit 
dependence on wi, so we need to know how it evolves. Unfortunately, the search for an evolution equation for wi is 
problematic as we discuss below. 



3. The We equation and higher order terms 

Equations (|43|l . (|48|l and (|5U|l depend on the orientation parameter wi. We interpret this as the degree to which the 
strings lie in the extra dimensions, that is, a hidden small-scale structure parameter. To make a fully closed system 
of equations we also need an evolution equation for wi. 

In analogy to our treatment of string velocities, we can try to obtain an evolution equation for wi by differentiating 
its definition H45|) with respect to r. Unfortunately this produces terms of the form 

I ((x'^Xi^) - (x'^i^)) , I ((x'^)' - (x'^)) , I ((x2)(l'Va^> - (i^l'VO) , etc. (58) 

These purely statistical terms are higher order and cannot be easily determined. We found terms of similar kind in 
the a:- velocity equation (|46|l . but in that case their only effect was to 'renormalise' the coefficients of other terms in 
the equation. Here, however, there are no such terms to be renormalised and the small differences between these 
undetermined statistical terms contribute at leading order. 

However, there are special cases in which an evolution equation for wi is not needed. For example, the evolution of 
a string network in a higher dimensional FRW universe is described by equations 1)41(1 and H57|l . which do not depend 
on wg. In particular, isotropy suggests that string segments would have equal probability of moving in any direction 
and one would expect to be a constant, namely the square root of the ratio of the number of extra dimensions, 
D — 3, over the total number of spatial dimensions D. 

For the case of small wi we note that since it appears in our equations only through factors of (1 — w|) and (2 — wf) 
modifying various coefficients, its evolution in time can be ignored as long as it stays small at all times. This is 
conceivable in cases of FRW universes with static extra dimensions, where the Hubble expansion of the three FRW 
dimensions will tend to reduce the value of wg. 

Finally, in cases that none of the above are applicable, it may be possible to find an ansatz for wg as a function 
of the different velocity components, based on physical arguments. We shall consider all three possibilities in the 
following discussion. 



General Features of Solutions 



(i) Isotropic case: We begin by considering the isotropic case b{t) =a{t) corresponding to a (_D + l)-dimensional FRW 
universe. The energy density evolution is described by equation H41|) . Comparing this to the corresponding equation 



in 3+1 dimensions H24|l we see that the effect of extra FRW dimensions is to modify two of the terms in (|24|l by a 
factor of 2/{D — 1). This reflects the fact that in three spatial dimensions, the density of the network is inversely 
proportional to the square of the correlation length L (equation H17|ll. whereas in I? > 3 dimensions it is inversely 
proportional to the (Z? — 1)'^ power of L (equation H35fl '). The energy density in the higher dimensional case decreases 
faster, as there are more expanding dimensions (also see equation H36() ). Of course there is a second, more dramatic, 
effect related to string interactions. As noted in 1111 B H the loop production parameter cd in the higher dimensional 
case is much smaller, since it is harder for two strings to find each other and intercommute in > 4. In fact, cd 
is proportional to (7^)"'^"^^ and as a result equation (|41|l does not have a scaling solution of the form 7 — const. 
In particular, the loop production term becomes smaller as time increases and 7 will always decrease, so that strings 
will dominate the energy density of the universe. This can be seen in Fig. ^were equations 141|l and H57(l with cd 
given by H4()|l have been solved numerically for Z? = 4. 

The asymptotic solution for L in a ZJ-dimensional isotropic model can be easily found to be 



Contrast this EDVOS model result with simple conformal stretching of the string network L oc a, obtained by assuming 
w = 0. In the EDVOS model if we start with v ^ 1 wc still find L oc a, but we also discover that v increases according 
to V oc t^^^ . Thus, conformal stretching is only a transient solution, which will be followed by one with non-negligible 
V. The corresponding scaling result is 

L cx t^^'^D-i oc a°^o-i , V = const . (60) 

In flat space cannot exceed 1/2 but in an expanding space this upper value is somewhat reduced. Since the horizon 
is expanding linearly in time, the correlation length quickly falls behind and the cosmological expansion becomes 
irrelevant so that 1/2. Equation (|60|l then implies (|59|l asymptotically. 

Note that (jHOJ means that the string energy density scales as oc L ^' oc a , which decays slower 

than radiation prad oc a^^^^^-*. Therefore, if this regime persists for long enough, the strings will eventually dominate 
the universe. The Friedmann equation for a string dominated D + I FRW universe yields 

aoct^/f^+^-'-i) (61) 

Substituting in H6()|l we finally obtain 

L cx . (62) 

The above reproduces and generalises the recent findings of Ref. for D = 3. In particular, we recover non- 

intercommuting strings in the linear scaling solution L (x t (i.e. setting P = and D = 3). 

In a situation were the extra dimensions are compactified, but expanding isotropically with scalefactor a(t), the 
solution H59|l implies that the correlation length will eventually catch up with the size of the extra dimension Ri and 
so the string network will become effectively three-dimensional with cd — c {6/Re) ^ where L > Ri (x a. Unlike the 
case of static compact dimensions, which we will discuss, this change is insufficient to prevent string domination over 
ordinary matter and radiation. 

(ii) Anisotropic case: We now turn to the case b(t) = 1. The relevant equations are (|43|l and (|48I50() with and 
ki approximately given by H54I55|I . To demonstrate the general features of solutions, we can keep Wi as a constant 
parameter (in the range < wi < 1), whose only effect is to modify the coefficients of damping terms in H43(l . H48|l 
and (|50|l . We will later consider a simplified version of the EDVOS model, where is not constant, but is given by 
a physically motivated ansatz, which is a function of vg and Vx- 

If the string network at formation is Brownian in D dimensions, then the curvature vector u explores the a;-space 
as well as the £-space and k^, kg are comparable. Equations H48() for Vx and H5()|l for vi have comparable source terms 
but ve has a much weaker damping term. Therefore, we expect Vi to become much greater than Vx and eventually to 
drive Vx to zero, because of the constraint Vx'^ + vi' < 1/2. We conclude that expansion of the x-space together with 
the fact that the ^-space is not expanding will halt string motion in the x-space. This can be verified numerically (Fig. 
EJ by solving (|48|l , (jSOfl with the approximations ^ ~ ^ and ^ ~ ^ in H54I55|I . These are valid approximations in 
both the non-relativistic and relativistic limits, as long as the curvature vector explores all D spatial dimensions. Thus 
the strings in this case will soon become non-interacting and, if this regime were to last long enough, dominate the 
universe irreversibly. However, there are two important caveats in this case, namely that anisotropic expansion would 
soon distort the simple Brownian network structure and that the intercommuting probability for a string network 
in D spatial dimensions would be small and decreasing (as discussed above). For compact extra dimensions, this 
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FIG. 1: Evolution of 7 in a log-log plot for a string network evolving in a (4 + l)-dimensional FRW universe. The capture 
radius 5 of the strings has been assumed to be equal to the initial correlation length and the exponent parameter /3 has been 
set to 1/2. There is no scaling solution and at late times log{'y) has a constant slope of —1/3, in agreement with our asymptotic 
solution l|59^ . 
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FIG. 2: Velocity evolution for a Brownian network in D spatial dimensions, three of which are expanding. The expansion 
strongly redshifts the velocities in the expanding dimensions {vx), resulting in domination of the extra dimensional velocities 
(ve). The strings will stop moving in the expanding dimensions. 



regime will end when the correlation length becomes larger than the compact dimension L > R(, after which point 
the evolution becomes effectively three-dimensional (we will return to this point later). 

The situation is different if the curvature vector effectively lies in the x-space, as may be the case in branc inflation 
(see next section) where the formation of the string network is constrained on an effectively (3 + l)-dimensional FRW 
slice. In this case the dot product u - 1 is negligible and thus is much smaller than k^- Equation H5U|) has essentially 
no source term. We expect to see Vx > vi with V£ given by its initial condition, weakly damped according to equation 
(I50II . This was verified numerically (Fig. ^ by setting ^ ~ as before, and ^ <C 1. The effect of increasing ^ 
was also considered. It was found that there is a critical value of ^ ~ 0.15 for which the scaling values of Vx and ve 
become equal. Above this critical value vi eventually dominates and we return to the previous regime, with ki and 
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FIG. 3: Velocity evolution for a Brownian network formed on a FRW slice. The curvature vector lies in the space of the three 
expanding dimensions, so there is no source for the extra dimensional velocities vi. They are given by their initial condition, 
redshifted weakly by the expansion. The three-dimensional velocity will not vanish but it can be significantly smaller than 
in the purely three-dimensional case i; ~ 0.7. 



kx comparable. In particular, for — ~ 0.5, soon reaches a 'relativistic' value, driving to zero. It may seem 
that if the strings are free to explore the extra dimensions after formation we may run into the same problems we 
had before, namely a practically zero intercommuting probability and the loss of the Brownian structure of the string 
network. However, the extra dimensions can be small and compact, in which case the probability of intercommuting 
can still be appreciable and the network is to a good approximation Brownian in three dimensions. 

Finally we consider the effect of varying the parameter wi. Wc find that 7 is insensitive to changes in wf between 
and 0.2 but is reduced by approximately 20% for wj — 0.5. On the other hand vi depends on wi more strongly. In 
particular, changing wf from to 0.01 increases the scaling value of vi by approximately 0.5% but a further increase of 
wl to 0.1 increases typically by 30% (this does not strongly affect Vx^ which is approximately equal to \Jl/2 — vi^). 
As a result, the effective three-dimensional behaviour of the evolving network is insensitive to the value of at least 
when it is smaller than 0.1 or so. There is a critical value of about 0.2 above which the sign of the damping term 
in equation (|50|l becomes positive and vi begins to increase. This apparently unphysical effect signifies the fact that 
the constant wi approximation breaks down for relatively large values of wi. In the following we discuss a simplified 
version of the extra-dimensional VOS model with wi given by a function of vi and Vx, where this unphysical behaviour 
does not occur. 



5. A Simplified Model 



We can obtain a simplified version of equations H43|l . (|48|l and (|50f) by expressing wi as a function of vi and Vx- 
String gradients in one direction are produced by velocities in the same direction so one would expect the network to 
achieve a sort of equipartition between velocities and gradients in the extra dimensions. This motivates the ansatz 



which mathematically corresponds to 



,.2 _ f^i! 

„2 



wi = ^ (63) 
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(64) 



With this substitution, equations H43|l and (|48I50() become 




(65) 



(66) 



(67) 



These three equations describe the macroscopic evolution of a Brownian network of strings, produced in a three- 
dimensional FRW slice, which are then free to explore the extra dimensions. In order for the Brownian structure to 
be maintained (so that equation (|65|l is valid) the extra dimensions need to be compactified at a size smaller than the 
initial correlation length at the time of formation. 

Note that the unphysical behaviour leading to the change of sign of the damping terms in the evolution equation 
(which is a result of the fact that the constant wi approximation is not valid for large wi) is now absent from equation 
(|67|l . The only way the damping terms can change sign is if exceeds 1/2, which is not allowed by the constraint 
(|56|l . Similarly the damping terms in equation (|66|l cannot change sign either. 

We can now solve (I65I67|I together with (|54I55|I for and k^. Since we are considering the case where the structure 
of the network at formation is essentially three-dimensional, the curvature vector lies mainly in the x-space, so that 
ki ^ kx (equivalently ^ ^ 1 in equation (|55|l ). The solutions are qualitatively the same as the corresponding 

constant wi solutions mentioned above. This is because wf = vt^ jv^ is small (and decreasing with time) so its effect 
in equations (|65I67|) is only to slightly modify coefficients of order unity. We will discuss these solutions in more detail 
in the next section where we will apply the model to branc inflation. 



Recently there has been much activity in trying to derive models of cosmological inflation from string theory 
[2l|2l|2alMIMl3l'33,l35l,l3lls3,Sl2a0 which there are a large number of scalar fields that can potentially 
serve as inflatons. One of the most interesting scenarios (especially from the point of view of cosmic string production) 
is brane inflation [27L 12^ l29l Isol Isij , where the role of the inflaton field is played by the distance between two branes 
or a brane and an anti-brane. The branes are initially relatively displaced and move towards each other due to an 
attractive potential arising from the exchange of bulk NS-NS and R-R modes. As they come closer together, open 
string modes stretching between them start contributing more strongly to the infiaton potential. At a critical distance 
of order the inverse superstring scale Mg"^ such modes become tachyonic and inflation ends in a hybrid-inflation- type 
exit. 

Brane-stacks carry Chan-Paton gauge groups and hence the tachyonic instability appearing at brane collision 
corresponds to a symmetry breaking process. This allows the formation of topological defects as lower dimensional 
branes if the vacuum manifold is non-trivial (SJ. For two brane-stacks of N branes each, the vacuum manifold is 
isomorphic to U{N), whose only non-trivial homotopy groups are the odd ones, n2k-i- Thus the topologically allowed 
defects have even codimension 2k. 

We need to make sure that cosmologically dangerous defects like monopoles or domain walls are not produced. This 
is fortunately the case (but see ^ can be seen by the following argument given in Ref. In order to describe 
our (3-|-l)-dimensional universe the branes must have three infinite spatial dimensions. They can either be D3-branes 
or Dp-branes with (p— 3) dimensions compactified at a size smaller than the horizon Thus the Kibble mechanism 
can only take place in the three infinite dimensions so the codimension of the defects must lie in these dimensions 
(hence it can only be 1,2 or 3). Since the codimension is even, the defects that are produced have codimension 2, that 
is they are D{p — 2)-branes wrapping the same compact space as the original Dp-branes. These objects will appear 
as cosmic strings to a three-dimensional observer. 

Alternative string production mechanisms (bulk preheating, resonant string formation) have been discussed in Q, 
resulting in both D-string and F-string networks, or even an interacting (F-D)-string network (also see ^3)- String 
stability was studied in 44]. Various possible cosmologically relevant metastable strings were found in many models, 
including {p, q) strings that is, bound states oip F-strings and q D-strings. For string stability also see |45l |. 

Given these possibilities, cosmic string production can be considered as a generic prediction of brane inflation 
(though models with no cosmic strings have been constructed I* is then natural to ask what the evolution of 

the string network would be in these models and whether it is distinctly different from the standard cosmic string 
evolution, thus providing a potential observational window into superstring physics. We consider this question below. 



IV. APPLICATION TO BRANE INFLATION 



A. Basics 



We have seen that in the context of brane inflation, D-strings are in general D(p — 2)-branes wrapping the same 
(p — 3)-dimensional compactified space as the initial Dp-branes. Their production is localized on the plane of brane 
collision, but later they are free to explore the (9 — p) compact dimensions transverse to the brane. 

For two strings to collide their worldsheets must intersect. In the 3+1 case the number of spacetime dimensions 
equals the sum of the dimensions of the two worldsheets and hence strings will generally collide . For a higher 
dimensional spacetime, however, this is no longer the case so strings will generally miss. One then expects to be able 
to model cosmic string evolution in the context of brane inflation by simply introducing an intercommuting probability 
P in the usual (3 + l)-dimensional evolution equations. 

This was done in I3l| where an intercommuting probability P (A/Ao in the original paper) was introduced in equation 
H19|l . The new scaling solution is 7 ~ P (cf 7 ~ 1 in three dimensions), corresponding to an enhancement in the 
scaling energy density of strings by a factor of P^^. This can be orders of magnitude different than the usual 3D case 

i 

There are at least two effects which could significantly alter the above result. The first is related to the small-scale 
structure of the strings. If they are wiggly, they may have more than one opportunities to intercommute during each 
crossing time. One expects that a probability of 1 /3 or so should have no impact in the scaling of 7, and thus, contrary 
to what is sometimes suggested in the literature, strings with such large P could not be observationally distinguished 
from standard field theory strings. On the other hand, for P <C 1 the scaling value of 7 should depend more strongly on 
P. The dependence of 7 on P, taking into account small-scale wiggles can only be modelled numerically. Simulations 
in flat space suggest a scaling 7 ~ \/P in the range 0.05 < P < 0.3 0, IH, To deal with this uncertainty 

we will introduce an effective intercommuting probability Peg = /(P) as a multiplicative parameter, modifying the 
phenomenological loop production term of equation H65(l . A numerical study of the functional dependence of P^s on 
P in an expanding universe will be presented in [5(^ . 

The second possibility arises from the observation that the velocities in the three infinite dimensions and in the 
compact ones, Vx and vi respectively, satisfy the constraint v^"^ + vi^ < 1/2. Together with the fact that Vi is very 
weakly damped compared to Vx fsection lIIIB4|l . this will slow down the motion of the strings in the three infinite 
dimensions. It is even possible that string motion in these dimensions can completely stop, in which case the strings 
will no longer be able to intercommute. We study this effect in more detail below. 



B. Applying the EDVOS Model 

As we have seen, the spatial structure of the string network at formation is essentially three-dimensional, but the 
strings evolve in 9 — p -I- 4 spacetime dimensions with 9 — p of them compactified. The compactification radius is 
assumed to be stabilised at a size smaller than the horizon [23| ■ Thus we can model the evolution of the string density 
using an effective three-dimensional description, where we introduce an intercommuting probability P to account for 
the fact that strings can 'miss' when they evolve in higher than 4 spacetime dimensions. As an evolution equation for 
7 we therefore use (I43II and replace the loop production term cvx^^^ by Pcgcvxj~^ , where an effective intercommuting 
probability Pcs = f{P) has been used to implicitly take into account the effect of string wiggliness on the scaling 
dependence on P. We also need to include both Vx and , the velocities in the (3 + l)-dimensional FRW x-space and 
the (9 — p)-dimensional compact £-space respectively, which we evolve using equations (|48ll . (|50|l . The parameters kx 
and kt are given by I|54I55|I for L» = 9 - p + 3 with ^ ~ ^ and ^ < 1 (see section Im B 4ll . 

The formation of the network takes place on the plane of the colliding branes, which has a finite thickness. One may 
worry that this fuzziness will translate to an uncertainty in the position of the strings in the extra dimensions, which 
may give rise to a significant structure in the £-space just after the strings are formed. This could spoil the assumption 
that the initial structure of the network can be considered three-dimensional, in which case equation (|43|l would not 
be valid. Fortunately this thickness is of the order of the inverse superstring scale , which is approximately 10"* 
times smaller than the correlation length at formation (see below). The network is to a very good approximation 
three-dimensional even at the time of formation. This justifies the assumption k^ <C kx or equivalently — ^ 1 at 
early times. 

The correlation length at formation, estimated by studying the tachyon potential, is of the same order of magnitude 
as the expected horizon size at that time 

La - H-^ - 10^ , (68) 
where the superstring scale is set by CMB data to a low GUT scale 



Ms - lO^"^ GeV. 



(69) 



0.6- 



0.5- 



vx, vl 



0.7 i 



0.7 H 



vx, vl 



FIG. 4: Velocity evolution in arbitrary time units with initial condition for the total velocity vo = 0.1, 1/y/TO, l/\/3 and l/\/2. 
We have assumed equipartition of kinetic energies ai t = and that the number of dimensions transverse to the brane is 2. 
For small veo the three dimensional velocities scale to a value close to the purely three-dimensional result v ~ 0.7. However for 
Vio > 0.4 the value of Vx can be significantly less than 0.6. 



After the strings are formed, the correlation length grows because of the Hubble expansion. On the other hand the 
dimensions transverse to the brane are compactified at a size a few times the inverse superstring scale, which is some 
10'^ times smaller than the correlation length. Thus, if the initial long string network is Brownian, it will remain so 
(as in the usual three-dimensional case) and equation H43|l will be valid at all times. 

\ simplify the equations we can use the ansatz voj = ^ (section IlIIB 5|l . Equations H43|l . H48|) and (|50|) are then 
replaced by 
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where we have explicitly written the loop production parameter in terms of the effective intercommuting probability 
Peff, that is, we have set 



C PoffC 



(73) 
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FIG. 5: Velocity evolution in arbitrary time units for Vx <^ vi &t t — 0. Here, Vx rapidly increases and reaches the same 
equilibrium value as in the case Vx ^ ve Si.t t — 0. This shows that the assumption of equipartition of kinetic energies at string 
formation is not necessary. The important initial condition is ve. 
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FIG. 6: Velocity evolution in arbitrary time units for vi — 0.95/\/2 at t 
about 0.2. Such a low velocity also affects the string density (Fig. [7|l. 



: 0. Because ve is so high, Vx scales to a low value of 



Below we discuss numerical solutions of equations (|7UI72ll with different initial conditions for Vx, vi as well as 
different values for the effective intercommuting probability Pcs in the range between 10^"^ and 1. For an illustration 
we first consider the case Poff — 0.1. As in section UlI B 41 we have that since ^ is small, there is essentially no source 
term for vi in (|72|l and thus vi is just given by its initial condition slowly damped by the expansion. On the other 
hand Vx is sourced by the string curvature R ^ L and although it is more strongly damped it dominates ve ■ There is 
a critical value of — ~ 0.15 above which the curvature (source) term for vi becomes large enough for ve to dominate. 
This value is too large to be relevant in the context of brane inflation: it corresponds to a situation where more 
than one tenth of the velocity developed during the last correlation time is in the ^-dimensions, which would require 
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FIG. 7: The effect of slowing down of strings on 7. On the left is the evolution of 7 for vi ~ 0.3 at formation. On the right is 
the corresponding evolution for vi — 0.95/v'2, in which Vx is considerably smaller. 



the curvature vector to have a significant component in the ^-space. As explained above the correlation length (the 
typical radius of curvature of strings) at formation is of order 10'^ times the size of the extra dimensions and further 
grows with the expansion, so — is expected to be much less than this critical value. Thus, string propagation in the 
x-space will not stop, but it can be significantly slowed down if the strings are created with enough momentum in the 
extra dimensions. 

Indeed we expect that some of the energy associated with the brane collision will be converted to kinetic energy 
of string segments in the dimensions transverse to the brane, so the slowing down of string motion in the infinite 
dimensions could be a significant effect. In Fig. ^we plot and vg as functions of time, assuming equipartition of 
kinetic energies at formation. This assumption is not necessary, and even if Vx ^ f£ initially, Vx will increase very 
fast (Fig. [SJ because of the curvature source term in equation (|71|1 . As long as the constraint < 1/2 is satisfied, 
the important initial condition is vi at the time of string formation. From Fig. ^we see that for initial values of vi 
smaller than 0.4 the scaling value of Vx is less than a few percent different from the purely three-dimensional result 
V ~ 0.7. However for vi initially greater than this value, the slowing down of strings in the three infinite dimensions 
is significant. In particular if vi has a value close to the maximum allowed {vi ~ 1/^/2) then Vx will approach a very 
small value Vx — ^1/2 — v^^ (Fig. ^ . Such a dramatic reduction of the three-dimensional speed of the strings also 
has a significant effect on the scaling value of 7 (Fig. ^ because the loop production term is proportional to Vx'- 
slowly moving strings will intercommute less often so the final density of the network will be higher, corresponding 
to a smaller 7. 

It is interesting to compare our results with the discussion of Jones, Stoica and Tye As a first approximation 
they used a simple one-scale model and assumed a relativistic, constant speed of strings. They found that the scaling 
value of the energy density of the network is times greater than in the purely three-dimensional case. The leading 
correction to this result comes from allowing the velocity of strings to be a variable, that is using a VOS model rather 
than a simple one-scale model. This is equivalent to setting = in our extra dimensional EDVOS model. The 
effect of this variable velocity is to reduce the density by a factor of 10 or so. We then include velocities in the extra 
dimensions vi ^ 0. If < 0.4 there is no additional observable effect. \ivi > 0.4 but considerably smaller than 1/V2 
then these extra dimensional velocities have no significant effect on the energy density, though they may lead to a 
substantial reduction of the three-dimensional string velocities. Finally, if ~ l/v^ the strings will be moving very 
slowly in the three infinite dimensions so that the number of intercommutings will be further reduced. This results 
in a substantial increase of the energy density of the string network. 

Finally we study the effect of varying the intercommuting probability P. For networks of the same type (F or 
D-string networks), P is expected to be in the range 10^'^ < P < 1 (in particular lO^'^ < P < 1 for F-strings and 
0.1 < P < 1 for D-strings !5|). However, since the strings can develop significant small-scale structure and become 
wiggly, they may have more than one opportunity for reconnection in each crossing time. Our one-scale model does 
not fully take into account such wiggly effects, but we have done so implicitly by using the effective intercommuting 
probability PcS = f{P) > P, as we have discussed above. The effect of small-scale wiggles can only be accounted for 
numerically and the results of such studies in an expanding space will be presented in . For the present discussion 




2000 4000 6000 SOOO 10000 



0000 20000 30000 40000 50000 



20000 40000 60000 SOOOO 100000 



FIG. 8: The dependence of 7 on the effective intercommuting probabihty Pcff- From top left to bottom right Pcff = 
1,0.5,0.3,0.1,0.01,10"^. As Pcff tends to zero the approximation 7 ~ Pcff becomes better but the time required to reach 
scaling increases. 



we take Pcff in the range 10^'^ < P < 1 and study how the model is affected by changes in Pcff- In Fig. |Slwe show 
the behaviour of 7 for different values of Pcff in that range assuming a moderate Vi ~ 0.36 at formation. Reducing 
Poft leads to less intercommutings and therefore greater string energy density, corresponding to a smaller 7. 



V. SUMMARY AND DISCUSSION 



We have presented an extension of the VOS model that can be used to study the macroscopic evolution of a 
Brownian network of Nambu-Goto cosmic strings in cosmological spacetimes with extra dimensions (EDVOS model). 
In order for the Brownian structure of the network to be preserved by time evolution either the spacetime must be 
isotropic or, in the case of anisotropic expansion, the network must be formed on an isotropic slice of the spacetime. 
If the strings are not confined on this slice after their formation then the extra dimensions must be compactified at 
a size smaller than the correlation length, as is the case in models of brane inffation. The evolution has then an 
effective, three-dimensional description, in which the effect of possible velocities in the extra dimensions can be taken 
into account. This is a significant factor because extra-dimensional velocities will act to slow down string motion in 
the infinite dimensions, reducing the number of string intercommutings while also changing the strings' effective 3D 
energy per unit length. 

First, we applied the model to the case of an isotropic [D + l)-dimcnsional FRW universe and found the generic 
behaviour (for D > 3) in which L (x a^/^^~^^ (in contrast to the naive conformal stretching obtained without the 
EDVOS model L (x a). Obviously, these strings do not scale because they find it increasingly difficult to find each 
other and intercommute, and so they would quickly dominate the energy density of such a universe. However, even 
if the additional dimensions were compact, this evolution might pertain at early times, as long as the expansion 
remained isotropic. In this case, the network correlation length L would continue to expand until it inevitably catches 
up with Ri oc a, the scale of the expanding compact dimension. Thereafter, the network would become effectively 
three-dimensional and, if the compact dimensions were then stabilised, the evolution would turn over towards the 
alternative FRW 3-brane scenario we discussed earlier. We point out therefore, in the case of static compact extra 
dimensions that, even if we begin with an initial correlation length smaller than the compact dimension Lq < Rg, the 
correlation length will grow until L > Ri, at which point the evolution becomes effectively three-dimensional. Hence, 
we conclude that the usual assumption Lq 3> Ri is not actually necessary for achieving asymptotically 3D scaling 
evolution. 



The most pertinent application of our model is to a situation where the string network is formed on a FRW 3-brane 
with the dimensions transverse to it compact and small (we have only considered flat extra dimensions). We have 
allowed the strings to be able to explore the bulk after formation, as is the case in models of brane inflation. We 
find that the density of the network in our VOS model (after scaling has been achieved) can be up to a factor of 10 
smaller than previous estimates Q. This correction comes from quantitatively accounting for the role of the string 
3D velocity, allowing it to be a variable rather than a constant as in the simplest one-scale models. We also found 
additional effects arising from the fact that the strings can also move in the extra dimensions. Standard field theory 
strings can be shown numerically to approach an average velocity of about 1/a/2 (on lengthscales A <C H~^). This is 
also the case in our extra-dimensional VOS model but now a significant amount of this velocity can be trapped in the 
extra dimensions, with the result that the observable 3D velocity is reduced. Since only the three FRW dimensions 
are expanding while the size of the extra dimensions is stabilised, there will only be a very weak redshifting of the 
extra-dimensional velocities. Thus, if the strings are created with significant velocities in the extra dimensions, these 
could survive for a long time and result in a slowing down of string motion in the three large dimensions. 

One could argue that after reheating, the strings would enter a damping regime where velocities would slow down 
due to the damping force from a high radiation background density. This could damp velocities in the extra dimensions 
more strongly than Hubble expansion, in particular eliminating v^. However, we note that reheating only takes place 
on the brane (the bulk remains cold and empty j27j ) and it is only the small portion of the string network intersecting 
the brane which feels this damping. As the strings move, they briefiy encounter the brane and pass through it, feeling 
a damping force, but at any time most of the length of the network is in the bulk, where no friction is felt. There 
are only a finite number of such brief encounters before cosmological expansion on the brane cools it sufficiently for 
the friction force to become negligible. The strings also might be oriented such as to intersect and remain in contact 
with the 3-brane in the longer term, but this only implies that the string will be 'pinned' at one point. Taking into 
account these geometrical effects, it appears that the usual 3D frictional damping terms will be suppressed by factors 
of order the intercommuting probability P (for both vi and Vx)- The network in higher dimensions will approach its 
relativistic scaling regime much more quickly than its 3-f 1 FRW counterpart, though a more quantitative investigation 
is certainly warranted. 

Gravitational back-reaction can also be expected to reduce the extra-dimensional velocities. The effect of gravita- 
tional backreaction can be incorporated in the VOS model by including a term SFG/iw^ in the evolution equation for 
the correlation length L [23|) where for long strings F is a constant of order 10 j2^. Even if G/i has a value close to 
the upper limit 10~^, this term is small compared to the loop production term Pcv so it has little direct infiuence on 
the network density. However, its effect on the extra-dimensional velocity could be more dramatic since gravita- 
tional radiation tends to act to eliminate small-scale structure below a certain length-scale A. This lengthscale was 
previously thought to be A ~ TGfit, but closer investigation has suggested that it was in fact considerably smaller 
A ~ {TG^Yt with Tj — 3/2 in the radiation era and rj = 5/2 in the matter era The reanalysis demonstrated 

that only modes of comparable wavelength interact efficiently, a fact which may be relevant in the case with extra 
dimensions. Here, we could envisage the small-scale modes trapped in the extra dimensions essentially decoupling 
from the effective 3D evolution. Damping of the velocity then would be primarily through self-interactions between 
the trapped modes which tend toward slow power law, rather than exponential, suppression. Whatever the outcome 
of a closer examination of this issue in higher dimensions, it is clear that the extra dimensional velocity vi will have a 
long lifetime and will influence the network evolution over long timescales. Note however that since the string position 
in the compact dimensions is described by (worldsheet) scalar fields, one expects some stabilisation mechanism for 
these moduli to kick in at low energies. This would render these excitations massive, effectively localising the string 
in the extra dimensions 44] . Given the large hierarchy between the preferred inflationary and SUSY breaking scales 
in brane inflation models (GUT and TeV respectively), the extra-dimensional effects we have considered are likely to 
be relevant for many orders of magnitude in time. 

While our EDVOS model can characterise the effect of the velocity ve on the overall large-scale 3D network prop- 
erties, uncertainties remain as to its magnitude at network formation, the subsequent (weak) rate of damping, and 
whether it can be sourced in any way during the subsequent evolution. Given the brane collision out of which the 
string network forms, it seems plausible that the string network will be released from the brane with a significant 
velocity vi with which to traverse the extra bulk dimensions. If V£ is relativistic, then the average 3D velocity Vx 
will be smaller and the network density will be higher than previously expected. However, this presumes the ab- 
sence of any significant source terms for V£, for example, due to string reconnections or other dynamical effects. By 
including such source terms, we have demonstrated that there is, in principle, the pathological possibility that the 
Hubble-damped 3D velocity Vx becomes negligible because of continuous contributions to the undamped ve- In this 
case the string network would dominate over radiation in the early universe. But we have also shown that there 
is a threshold {ki 0.1^:^:) below which this does not happen. For brane inflation scenarios in which the string 
curvature is predominantly three-dimensional at formation, we expect this criteria to be fulfilled, with the EDVOS 
model predicting a scaling solution in which L = and 7 ^ Poff- We stress that the statement L o: P should not 



be taken as a prediction of our model. Instead, we have used an effective intercommuting probability Pcff, which is 
a phenomenological parameter of the model, to be determined by string network simulations or a deeper statistical 
analysis. The functional dependence of Poff = fiP) on the actual intercommuting probability P will be the subject 
of a forthcoming publication [H^ • 

Another interesting issue which we can use our EDVOS model to investigate is the evolution of closed strings 
(loops), and in particular the possibility that such loops could wrap around the compact dimensions 0. If the 
compact manifold admits non-trivial one-cycles, this would give rise to stable monopole-like objects (from the 3D 
point of view) which could dominate the universe or even provide a dark matter candidate. We denote these objects 
cycloops, that is, loops wrapping around non-trivial cycles. We leave the study of the formation, properties and 
consequences of cycloops for a different publication . 

The purpose of the present paper was not to exhaustively explore all possible avenues or to investigate specific 
brane inflation models, but rather to provide a broad picture of the key dynamical effects which will influence a 
cosmic string network emerging from a higher dimensional theory. We believe the generalised VOS model we have 
presented here is an important step in developing a more quantitative description of cosmic string evolution in brane 
inflation and other contexts. There are a host of further issues to explore. We have, for example, only considered the 
'abelian' case where there is only one particular type of string, which simply exchanges partners and reconnects when 
it meets another (albeit with a reduced probability for such an encounter). Similar methods can be applied to the 
richer structures created by F and D strings in more complex (true) string networks (see Ref. '3-2]) ^nd this deserves 
further investigation. 



APPENDIX A: APPROXIMATE FORMULAE FOR k, AND ki 

Deflning y= (x, 1/a) we split the velocity y into a 'curvature' component yc produced during the last correlation 
time and a 'left-over' component yp (coming from previous accelerations) by writing y = Yc + Yp- We interpret the 
flrst as the velocity induced on large-scales by the present correlation length of the string, whereas the second is the 
velocity remaining from previous correlation times, generally on small scales. These two components are uncorrelated 
so we have (yc • yp) = 0. We also have the gauge condition y ■ y' = and we assume that both components yc and yp 
separately satisfy this, i.e. we have yc • y' = and Yp ■ y' = 0. 

Hence, since the curvature vector u is normal to the tangent y' we have 

D-2 

Vi^Ai^- + J2 <^»y^ (Al) 

i=l 

where y,. , yp are the unit vectors in the direction of yc , yp respectively and y^ are unit vectors spanning the 
(D — 3)-dimensional subspace normal to yc , yp and y'. Dotting with y and taking the average along the string we 
have 

,^ I . 1 X \ / A f^ - B\yjA 



(u . y) ^ (A|yc| - S|yp|) = {^A\y^\ ~ ) Vy' {<fi + y^)'^'"^ - ( Vr ^ ) ■ (A2) 

1/2 

Taking the modulus of ljAl|l gives 1 = a(^1 + ^ + J2iLi^ ^2") > which we use to substitute for A in ljA2|l . 
Remembering that few ~ (u • y) we find 

Av, (A3) 



,1/2/. „n_3c.^^^/' 



For small velocities {v <^ 1) we have that i; ~ t^c and the ratio ^ is much less than unity. However, as the velocity 
tends to relativistic values, the Vp contribution becomes more and more significant. Assuming that this relative 
contribution is proportional to some power of the total velocity, we set ^ = fv" with /, a constants (clearly a > 1). 

Similarly in the low velocity limit, the curvature vector u is parallel to yc so ^, ^ <C 1 but as v increases, B and Ci 
become comparable to A. Assuming a power law dependence we set — = gv^ ■ Equation ljA3p becomes 



A' A 

A ~ A 

1 - fgv°'+'^ 



(1 + /2w2a)l/2 (1 + (D - 2)521,2/3)1/2 



(A4) 



Comparison with the hehcoidal string solution in D=3 flat space gives a + /3 = 3 |2J|. The limit fc(l/\/2) = 
(which can be shown to hold analytically p^l, "271) gives fg = 8. Then, comparison with the three-dimensional case 
(|3()|l suggests f = g, a = (3. The approximate formula for k is therefore 

(A.) 



(1 + 8^6)1/2 (i + 8(i?-2)t;6)i/2 
For kx we form (u • x) and work as above to find 



1 _ B Vxp Va 

Vxc ^ A 1,^,2 ■ 



— . s 1/2 / " s 1/2 (^^^ 



where we have used that (y^ • x) = (^j^^Yc ■ = <^ |J;|Xc • = {'WJ / 

Assuming that has the same velocity dependence as ^ we obtain the following approximate formula for k^ 
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k ^^'^ 1 — 8u (A7) 

^~ (1 + 8w6)^/^(l + 8(i:'- 2)^6)1/2 ■ 



The corresponding equation for kg is 



-Vc (1 + 8v^)^^^ (1 + 8(13 - 2)i;6)^/^ ' 
Equations (|A5|I . HA7|I and HA8p are indeed the approximate formulae (|53I55|I given in section UlTI 



(A8) 
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